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ABSTRACT 


It  is  desirable  to  perform  a  study  of  the  elastic 
responses  of  submarine  hulls  and  internally  attached 
structures  to  underwater  shock  loading.  In  view  of  the 
expense  and  possible  involvement  of  non-shock  testable 
equipment  it  is  not  always  desirable  to  perform  actual  shock 
tests  utilizing  real  hulls  or  SSTV's  (Submerged  Shock  Test 
Vehicles).  With  the  advent  of  large  scale  computing  power, 
numerical  methods  now  exist  to  predict  equipment/hull 
responses.  The  ELSHOK,  (ELASTIC  SHOCK)  computer  code  is 
used  to  perform  a  parametric  study  on  a  submerged  shell  with 
an  internally  attached  substructure.  Of  particular  interest 
is  the  phenomenon  of  dynamic  amplification  or  lack  of  such 
in  the  response  of  this  particular  substructure  and  effects 
on  substructure-shell  interaction  due  to  varying 
substructure  mass  and  stiffness  subjected  to  underwater 
explosion  side-on  shock  loading. 
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I.  INTRODUCTION 

A.  BACKGROUND 

A  submarine  must  be  able  to  withstand  moderate  to  severe 
shock  loadings  that  will  result  from  underwater  explosions 
(UNDEX)  it  might  be  subjected  to  in  modern  warfare. 
Conventional  and  nuclear  weaponry  deliver  devastating  forces 
and  no  vessel  can  be  expected  to  survive  in  the  events  of 
near  or  direct  hits.  The  emphasis  on  design  and  testing  is 
in  the  area  of  survival  of.  mission  critical  equipment, 
machinery  and  weapons  systems  of  a  platform.  The  ability  of 
a  naval  vessel  to  carry  out  its  mission  after  being 
subjected  to  an  UNDEX  will  depend  on  the  survivability  of 
these  installations. 

The  current  specifications  for  building  ships  and 
submarines  contain  the  requirements  for  shock  loading  to  be 
met  by  the  builder  or  vendor  of  shipboard  installations 
and  equipment.  In  general,  all  critical  equipment  is 
required  to  pass  a  series  of  shock  tests  where  testing  is 
practical  as  outlined  in  MIL-S-901D  fPef.  11.  This  document 
specifies  the  shock  testing  requirements  for  shipboard 
machinery,  equipment  and  systems  which  must  resist  High 
Impact  (HI)  mechanical  shock.  The  purpose  of  these  tests  is 
to  determine  the  suitability  of  machinery,  equipment  and 


9 


systems  for  use  after  exposure  to  severe  shock  which  may  be 
incurred  in  wartime.  It  is  not  always  practical  to  test 
equipment  due  to  size  and  weight.  In  the  case  of  equipment 
design,  it  is  not  always  practical  to  evaluate  an 
installation  prior  to  final  design.  Mechanical  shock 
testing  of  individual  equipments  may  not  always  indicate  how 
equipment  will  perform  or  respond  when  actually  installed. 

An  analysis  for  a  given  design  is  specified  for  the 
situations  where  testing  cannot  be  conducted.  It  is  based 
on  limited  and  dated  testing  and  observations  which  are  used 
to  specify  a  shock  desigin  factor.  A  static  design  force  is 
arrived  at  for  various  weights  of  equipment  which  is  used  to 
evaluate  the  design  of  support  structures  of  installations. 
The  method  is  simple  to  apply  but  does  not  always  correlate 
well  with  actual  testing. 

Numerical  methods  have  been  developed  for  use  in 


submarine  shock  response  analysis.  The  computer  code  ELSHOK 
(Elastic  SHock)  developed  by  Weidlinger  Associates  for  the 
Office  of  Naval  Research  [Ref.  2]  is  used  to  conduct  this 
study . 

B.  PURPOSE 

ELSHOK  has  been  developed  through  a  series  of  controlled 

tests  and  has  been  used  to  predict  equipment  responses  prior 
to  low  level  explosive  tests  on  submarines.  Although  ELSHOK 

has  been  validated  by  extensive  testing  it  is  considered 
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prudent  to  conduct  additional  research  to  gain  insight  to 
the  sensitivity  of  the  code's  predictions  when  a  variety  of 
conditions  are  applied.  The  intent  of  this  investigation  is 
to  examine  the  change  in  the  equipment/hull  interactions 
under  these  conditions.  The  models  used  to  conduct  the 
parametric  study  are  similar  to  those  used  in  the  mid  1970's 
in  tests  conducted  by  the  Office  of  Naval  Research  (ONR) 

[Ref.  3].  The  internal  structure  used  was  specifically 
designed  to  reduce  the  effects  of  dynamic  amplification  for 
the  conditions  applied  in  these  tests.  The  same  model 
geometry  was  studied  here  by  varying  one  main  dimension  to 
tune  the  substructure  to  the  natural  frequencies  of  hull  to 
see  if  this  particular  internal  structure  configuration 
retained  this  design  feature  over  a  range  of  mass  and 
stiffness  variations. 


UNDEX  SHOCK  RESPONSE  USING  ELSHOK 


A.  ELSHOK  PRINCIPLES  OF  OPERATION 

The  ELSHOK  computer  code  explained  in  [Ref.  2]  is  a 
suite  of  computer  programs  assembled  and  developed  to 
calculate  the  submerged  shock  response  of  a  linearly 
elastic  structure  subjected  to  an  underwater  shock  wave. 
Specifically,  a  ring-stiffened  shell  of  revolution  of  finite 
length,  with  or  without  internal  structure,  is  considered  to 
be  immersed,  initially  at  rest,  in  an  infinite  acoustic 
fluid  and  to  be  excited  by  an  acoustic  shock  wave.  The 
shock  wave  originates  from  an  UNDEX  arbitrarily  located  from 
the  structure  and  is  assumed  to  propagate  as  a  spherical 
wave  through  the  medium.  This  assumption  allows  incident 
fluid  particle  velocities  impinging  on  the  wet  surface  of 
the  immersed  structure,  the  hull  of  a  submarine,  to  be 
determined  from  the  incident  fluid  pressures.  ELSHOK 
contains  two  possible  representations  of  incident  pressure. 
One  uses  a  known  pressur e-t ime  profile.  This  type  is  of  the 
kind  that  may  be  used  to  describe  a  tapered  charge  simu¬ 
lating  a  nuclear  explosion.  The  other  is  the  empirical 
decaying  exponential  put  forth  by  Cole  [Ref.  4]  for  spheri¬ 
cal  charges  of  conventional  explosives.  If  the  shell 
contains  internal  equipment,  a  substructuring  technique 


described  by  Ranlet  [Ref.  5]  treats  the  internal  equipment 
response  by  coupling  the  free-free  modes  of  the  empty  ring- 
stiffened  shell  and  the  fixed-base  modes  of  the  internal 
equipment  through  use  of  dynamic  boundary  conditions. 
Structure-fluid  interaction  is  calculated  using  the  Doubly 
Asymptotic  Approximation  (DAA)  presented  by  Geers  fRef.  6]. 
Normal  fluid  displacement  of  the  structure-fluid  interface 
is  expanded  in  a  series  of  surface  expansion  functions  which 
are  orthogonal  over  the  entire  wet  surface  of  the  submerged 
hull  or  shell  model.  The  surface  expansion  functions  lead 
to  matrices  in  which  the  elements  are  determined  by  matching 
the  limits  of  zero  and  infinite  frequencies  of  the  pressure- 
velocity  relationships.  Thus,  exact  solutions  are  obtained 
for  the  transient  problem  at  early  and  late  times  and  DAA 
results  in  a  smooth  transition  between  the  two  limits. 
Because  DAA  calculates  the  effects  of  the  fluid  in  the 
interaction  problem  through  the  series  of  surface  expansion 
functions  defined  only  at  the  wet  surface,  the  fluid  field 
is  in  effect  uncoupled  from  the  structural  field.  The  fluid 
effects  produced  by  DAA  are  simply  applied  as  additional 
loading  in  the  modal  component  structural  analysis  of  the 
shell . 

The  substructuring  technique  separately  determines  the 
vibration  modes  of  the  components  that  comprise  the  combined 
submerged  structure  of  shell  and  internal  substructures. 
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Equations  of  motion  of  the  entire  structure  are  obtained  by 
combining  interaction  forces  and  moments  and  enforcing 
compatibility  of  deformation  at  component  connection  points. 
This  precludes  calculations  involving  the  modes  and  natural 
frequencies  of  a  system  of  combined  components  and  also 
avoids  a  requirement  for  a  combined  system  stiffness  matrix. 
The  ELSHOK  code  also  handles  the  case  of  no  substructure  in 
the  shell  using  the  same  substructuring  method. 

ELSHOK  is  particularly  suited  to  the  task  of  a 
parametric  study  due  to  its  purely  numerical  nature  and 
modular  organization.  Various  internal  equipment  models  can 
be  combined  repeatedly  with  the  same  shell  model.  Inputs  of 
pressure-time  profiles  or  decaying  exponential  types  are 
mathematically  applied  to  these  numerical  models  and  result 
in  output  in  the  form  of  velocity-time  histories.  This 
transient  response  analysis  takes  into  account  interactions 
of  shell  and  substructure  and  eliminates  dependency  on  data 
derived  from  actual  explosion  testing. 

B.  ORGANIZATION  AND  IMPLEMENTATION  OF  ELSHOK 

The  implementation  of  ELSHOK  is  accomplished  using  a 
suite  of  seven  major  programs.  These  major  components  are: 

1)  B0S0R4  -  structural  analyzer  for  shell  [Ref.  7] 

2)  ACESNID  -  virtual  mass  processor 

3)  PIFLASH  -  shell-fluid  processor 


4)  SAP  IV  -  structural  analyzer  for  substructure  [Ref.  P] 


5)  PICRUST  -  substructure  processor 

6)  USLOB  -  time  integration  processor 

7)  PCJSLOB  -  plotting  processor 

Illustrated  in  Figure  1  are  the  general  relationships 
between  the  main  components  and  four  phases  into  which 
ELSHOK  is  separated.  In  any  analysis  using  ELSHOK  the  shell 
model  and  internal  structure  model  must  be  available  or 
constructed  by  the  user.  Development  of  the  models  used  in 
this  study  is  discussed  later.  Once  models  are  arrived  at 
the  analysis  can  be  carried  out  in  the  four  phases  by 
executing  each  of  the  seven  major  programs  sequentially  in 
the  order  listed  above. 

1.  Phase  I — Shell  and  Fluid  Analysis 

The  finite  element  difference  code  BOSOR4  is  used  to 
analyze  the  idealized  model  of  the  submerged  hull  or  shell. 
The  model  must  be  symmetric  about  its  longitudinal  axis  and 
it  may  contain  rings  and  bulkheads.  Circumferential  rings 
can  be  treated  as  discrete  stiffeners  or,  if  closely  spaced, 
they  may  be  represented  in  the  orthotropic  approximation 
(smeared)  [Ref.  9J.  An  entire  main  body  or  full  model  is 
used  to  capture  the  gross  effects  of  the  shell  response  when 
calculations  are  made  for  circumferential  wave  numbers  N  =  0 
and  1.  These  include  rigid  body  translation  and  whipping  of 


the  structure  (N  =  1)  and  the  torsional  modes  of  the  struc¬ 
ture  (N  =  0).  Performing  the  analysis  for  N  =  0  and  >  2 


(breathing  mode  N  *  0)  obtains  the  local  shell  response.  In 
this  analysis  wave  numbers  N  *  0,  1,  2,  and  3  are  used. 
Figure  2  shows  the  basic  nodal  patterns  for  N  *  1-4.  ELSHOK 
requires  separate  runs  of  BOSOR4  for  each  wave  N  used  in  the 
analysis.  Each  BOSOR4  run  determines  the  in-vacuo  free-free 
modes  and  corresponding  natural  frequencies  from  the  dynamic 
properties  and  geometry  supplied  in  the  model  of  the  shell. 
If  the  shell  is  symmetric  about  the  midpoint  of  its  length 
and  the  shock  wave  impinges  on  it  so  as  to  produce  symmetric 
reponse,  only  half  of  the  structure  need  be  modeled.  Such 
is  the  case  in  this  study  since  the  shock  wave  emanates  from 
an  origin  located  normal  to  the  shell  at  midships. 

ACESNID  (Accession  to  Inertia  and  Damping)  accom¬ 
plishes  the  second  step  of  Phase  I  by  calculating  the 
virtual  mass  array.  This  provides  the  late-time  contri¬ 
bution  of  the  Doubly  Assymptotic  Approximation  [Ref.  6]. 
Normal  d is-pl acements  corresponding  to  the  surface  expansion 
functions  are  applied  to  the  surface  of  a  cavity  in  the 
fluid  (infinite  medium)  having  the  same  shape  and  size  as 
the  wet  surface  of  the  shell  model  being  investigated.  The 
cavity  geometry  is  obtained  from  the  shell  geometry 
contained  in  the  BOSOR4  output  files.  Only  one  run  of 
ACESNID  need  be  performed  to  encompass  all  values  of  N  being 
used  and  all  values  of  N  being  used  in  B0S0R4  must  be 
indicated  in  ACESNID  for  compat ibl i ty . 


PIFLASH  (Prepare  InPut  for  Fluid  jnd  Shell) 
accomplishes  the  last  step  of  Phase  I  and  that  is  to 
reorganize  the  fluid  and  shell  data  files  from  B0S0R4  and 
ACESNID  into  a  fluid-shell  file  that  is  used  to  solve  the 
response  equations.  If  there  are  no  internal  substructures 
involved  in  the  analysis,  then  this  fluid-shell  file  is  the 
main  input  for  Phase  III  and  Phase  II  is  not  necessary. 

2 .  Phase  II — Structural  Analysis  of  Internal 

Substructures 

A  finite  element  model  for  each  internal 
substructure  must  be  developed.  The  SAP  IV  code  performs 
the  structural  modal  analysis  on  each  substructure  to  be 
included.  The  code  contains  a  variety  of  types  of  elements 
that  can  be  used  with  ELSHOK.  Concentrated  masses  attached 
to  the  substructure  may  also  be  represented  by  the  methods 
used  in  SAP  IV.  Each  substructure  requires  a  separate 
performance  of  calculations  in  which  the  fixed-base  mode 
shapes  and  natural  frequencies  for  each  are  determined. 

This  data  is  recorded  in  the  "substructure  mode  file"  along 
with  the  geometry,  stiffness  and  mass  information  of  the 
model  for  future  processing. 

PICRUST  (Prepare  l_nput  for  Calculating  Response  of 
Substructure)  rearranges  the  data  from  the  "substructure 
mode  file"  and  computes  additional  information  required  to 
solve  the  response  equations.  PICRUST  also  provides  for  the 
selection  of  which  modes  are  to  be  included  from  SAP  IV  and 


the  specifics  of  connectivity  between  the  shell  and 
substructure  used  in  the  response  calculations. 

Additionally,  influence  coefficients  corresponding  to  the 
base  motions  of  the  substructure,  constraint  modes  and  modal 
coefficients  for  calculation  of  interaction  forces  with 
moments  between  the  shell  and  substructure  are  determined 
[Ref.  9],  PICRUST  is  the  last  step  in  Phase  III  and  its 
results  are  stored  in  the  substructure  file. 

3.  Phase  III — Shock  Response  and  Time  Integration 


The  information  required  to  complete  the  fluid- 
structure  interaction  is  supplied  as  one  input  file. 
Therefore,  if  a  substructure  is  being  analyzed  in  the 
problem  the  shell-fluid  file  from  Phase  I  must  be  combined 
with  the  substructure  files  produced  in  Phase  II  and  become 
the  input  file  for  USLOB  (Underwater  Shock  Loading  of 
Bodies)  which  is  ELSHOK's  time  integration  processor.  For 
the  case  of  an  empty  shell  the  shell-fluid  file  is  the  input 
for  USLOB.  Integration  is  accomplished  using  a  modified 
Runge-Kutta  scheme.  The  time  step  and  shock  wave  loading  is 
specified  by  the  user  in  USLOB.  Care  must  be  taken  to  be 
sure  that  the  product  of  the  time  step  chosen  and  the 
highest  circular  structural  frequency  <*>  of  the  modes  being 
used  from  SAP  IV  is  less  than  unity.  Otherwise,  the 
numerical  integration  scheme  will  not  converge  without 


changing  the  time  step. 


The  cases  studied  in  this  work  used  both  methods 


available  to  represent  the  incident  fluid  pressure.  A 
pressure-time  history  profile  shown  in  Figure  3  was  supplied 
as  a  typical  tapered  charge  simulation.  Nine  discrete 
pressure-time  points  are  entered  as  input  to  describe  the 
incident  shock  wave  loading.  Tapered  charges  are  used  in 
UNDEX  testing  to  produce  the  same  type  of  loading  resulting 
from  an  underwter  nuclear  detonation.  Here  the  impulse  of 
the  shock  wave  is  delivered  to  the  target  over  a  longer 
period  of  time  compared  to  the  interval  over  which  the 
impact  of  a  conventional  explosive  charge  is  delivered. 

Both  are  assumed  to  emanate  from  an  origin  far  enough  away 
so  as  not  to  include  the  target  in  the  gas  bubble  effects. 
ELSHOK  is  used  for  elastic  response  analysis  only  and  hulls 
within  the  immediate  vicinity  of  the  gas  bubble  formed  in  an 
UNDEX  would  most  likely  have  responses  well  into  the  plastic 
response  range.  Figure  4  shows  the  pressure-time  profile 
used  to  represent  the  incident  shock  loading  of  a  spherical 
charge  of  HBX-1,  a  conventional  explosive,  of  a  weight  which 
would  deliver  the  same  impulse  over  the  same  period  of  time 
as  the  duration  of  the  tapered  charge  loading  and  detonating 
at  the  same  range  from  the  closest  point  of  the  target  or 
standoff.  This  profile  is  an  empirical  exponentially 


decaying  pressure-time  history  where  the  initial  pressure  is 
the  peak  incident  pressure  on  the  shell  due  to  the  charge 


m 


weight  and  standoff.  It  is  seen  from  Figure  4  that  the 


majority  of  impulse  of  the  HBX-1  blast  is  delivered  inside 


the  first  millisecond.  The  empirical  equation  over  which 


USLOB  performs  integration  to  obtain  the  fluid  particle 


velocities  is  of  the  form  [Ref.  4] 


Pr(R,t)  =  K1(W-l/j/R)  exp(-t/0Q) 


where , 


Pj  (R,t)  =  the  incident  pressure  on  the  hull  at  range  R  (ft) 
from  explosive  source  to  point  of  interest  (psi) 


=  multiplicative  constant  for  incident  pressure 


=  spatial  decay  constant  for  incident  pressure 


=  time  after  arrival  of  shock  wave  at  point  of 
interest  (msec) 


*  weight  of  spherical  charge  (lb.) 


=  K3W1/3  (W1//3/R)  4  ,  time  constant  of 
exponential  decay  (msec) 


=  multiplicative  constant  for  time  constant  of 
exponential  decay 


=  spatial  decay  constant  for  time  constant  of 
exponential  decay 


USLOB  creates  a  generalized  velocity  file  for  the 


shell-substructure  and  shell-substructure  transformation 


file  to  be  used  in  the  last  phase  of  ELSHOK  in  producing 


plotted  and  printed  velocity-time  profiles.  The  trans¬ 


formation  file  contains  information  required  to  describe  the 


motion  at  the  shell-substructure  connection  points. 
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PUSLOB  (Plots  for  USLOB)  processes  the  results  of 
the  previous  three  phases  so  they  may  be  plotted  or  printed 
as  veloci ty-time  histories.  Velocity  time-profiles  are 
requested  for  various  locations  on  the  shell  and 
substructure  by  specifying  a  call  out  point  in  the  input 
files.  The  call  out  points  are  the  nodes  that  are  used  to 
describe  the  models  of  the  shell  and  substructure  in  BOSOR4 
and  SAP  IV.  Selection  of  nodes  is  important  because 
velocities  can  only  be  obtained  at  those  locations.  The 
shell  velocities  may  be  represented  in  either  a  global  or 
local  coordinate  system.  The  global  is  represented  by  X,  Y, 
Z  that  locate  the  positive  sense  of  X  toward  the  bow  or 
forward  section  of  the  shell,  the  positive  sense  of  Y  is  to 
port  and  positive  Z  is  up.  The  local  system  is  described  by 
u-longitudinal  (fwd  pos)  v-c i rcumf erenti al ;  w-inward  normal 
from  the  shell.  The  substructure  is  always  described  by  the 
coordinate  system  specified  in  its  SAP  IV  finite  element 


model . 


zation  of  the  ELSHOK  Code  [Ref.  2] 


circumferential  nodal  pattern 


Figure  2.  Nodal  Patterns  for  Cylindrical  Shells 


MODELS  INVOLVED  IN  THE  ANALYSIS 


A.  MODELING  WITH  FINITE  ELEMENT/FINITE  DIFFERENCE  METHODS 

Finite  element/finite  difference  methods  have  become 

powerful  tools  in  engineering  with  the  advent  of  large  scale 
computing  power.  These  methods  are  of  great  use  in  modeling 
techniques.  Even  though  the  models  are  much  simplified  for 
use  in  these  techniques  the  results  arrived  at  can  be  quite 
accurate  or  at  least  adequate  for  the  problems  that  might  be 
impossible  to  solve  by  other  methods.  Accuracy  of  the 
results  is  closely  tied  to  how  well  the  model  represents  the 
actual  configuration,  so  care  must  be  taken  when 
constructing  the  model.  In  ELSHOK,  the  finite  difference 
code  BOSOR4  is  used  to  describe  the  dynamic  behavior  of 
underwater  hulls  by  representing  them  as  much  simplified 
models  of  ring-stiffened  shells  of  revolution.  Similarly, 
the  models  used  to  analyze  internal  substructures  are  also 
simplified  representations  made  with  the  finite  element  code 
SAP  IV. 

B.  SUBMERGED  SHELL  MODELING 

Figure  5  is  a  schematic  of  the  general  shell  structure 


and  substructure  model  depicting  the  notation  and  convention 
used  to  describe  the  shell  global  coordinate  frame  X,Y,Z, 
shell  local  coordinate  frame  u,v,w  and  the  substructure 


coordinate  frame  x,y,z.  The  z-axis  orientation  of  the  local 
substructure  coordinate  system  is  related  to  the  Z-axis  of 
the  shell  by  the  angle  aQ  where  <*  represents  the  substruc¬ 
ture.  Points  are  located  along  the  meridian  of  the  shell 
using  the  variable  s.  The  angle  0  is  used  to  locate 
circumferentially  the  desired  local  shell  displacements. 

One  shell  model  is  used  in  all  calculations  of  this  study. 
Once  the  symmetric  full  model  was  constructed  no  changes 
were  required  to  conduct  the  investigation  of  varied 
substructure  sizes  and  types  of  UNDEX. 

The  shell  or  hull  that  has  been  modeled  for  analysis  by 
ELSHOK  in  this  study  is  similar  to  an  actual  structure 
configured  for  testing  by  the  Office  of  Naval  Research 
[Ref.  3].  Figure  6  is  a  schematic  of  the  representation  of 
the  full  shell  model  used.  Though  its  size  is  much  smaller 
it  is  also  of  similar  geometry  to  that  of  SSTV's  (Submerged 
Shock  Test  Vehicles)  that  are  used  to  investigate  the 
underwater  shock  effects  on  submarine  installed  equipment. 
The  model  is  an  idealized  submarine  hull  in  that  both  can  be 
considered  as  free-free  ring-stiffened  cylinders  with  end 
closures.  The  shell  in  this  study  is  in  fact  a  cylinder  of 
revolution  made  of  high  strength  steel  stiffened  by  six 
large  steel  rings  plus  forty  small  stringers  also  of  steel. 
The  end  plates  are  aluminum.  An  internal  substructure,  also 
similar  to  one  involved  in  the  testing  by  ONR  is  mounted  at 


26 


midships  on  the  port  side  of  the  hull.  In  B0S0R4  the  model 
is  divided  into  segments  of  consistent  properties  contained 
in  each  cross  section.  When  the  structure  is  symmetric 
about  the  midpoint  of  its  longitudinal  axis  only  one  half  of 
the  model  need  be  configured.  This  is  accomplished  by 
specifying  the  appropriate  boundary  condition  at  the  point 
of  symmetry.  In  this  problem  symmetry  is  employed  allowing 
the  shell  model  to  be  divided  into  two  segments.  Segment 
one  is  the  aft  aluminum  end  plate  and  has  nine  nodes 
specified  on  its  surface.  Node  1/01,  node  one  on  segment 
one,  is  the  center  of  the  end  plate.  The  other  nodes  move 
as  concentric  circles  to  the  outer  circumference.  Segment 
two  is  the  cylinder  from  the  end  plate  to  the  midships 
point.  Twenty-seven  nodes  are  specified  for  segment  two 
that  are  actually  bands  around  its  circumference.  Node 
numbers  increase  from  aft  to  forward  such  that  node  2/27  is 
the  midpoint  on  the  cylinder. 

Concentrated  masses  and  substructures  internal  to  the 
shell  must  be  modeled  separately  if  they  contribute  signi¬ 
ficantly  to  the  stiffening  of  the  shell.  Therefore,  the 
large  rings  have  been  modeled  as  discrete  stiffeners  and  the 
internal  substructure  is  modeled  separately  using  SAP  IV. 

If  there  were  no  interest  in  the  substructure  response  it 
could  be  modeled  as  a  concentrated  mass  attached  to  the 
shell.  There  is  only  one  substructure , a  doubly  overhung 


cantilever  beam,  used  in  this  study  and  there  are  no 
concentrated  masses.  The  smaller  stringers  do  not 
contribute  as  significantly  to  the  response  of  the  shell  as 
do  the  large  rings  so  they  are  smeared  onto  the  shell  to 
account  for  their  mass  [Ref.  9]. 

The  cylindrical  section  of  the  shell  is  a  surface  of 
revolution  having  consistent  properties  throughout  its 
entire  length.  This  allows  information  regarding  Young's 
modulus,  Poisson's  ratio,  mass  density,  thermal  expansion 
coefficients,  plating  thickness  and  whether  or  not  rings  are 
smeared  to  be  specified  at  only  one  point  for  the  segment  in 
BOSOR4.  Similar  information  is  specified  for  segment  one, 
the  end  plate.  Tables  I  and  II  are  listings  of  the  natural 
frequencies  and  wave  numbers  for  each  mode  calculated. 
Twenty-two  N  =  0  breathing  modes,  thirty  N  =  1  rigid  body 
and  whipping  modes  and  thirty  N  =  2,3  modes  were  retained 
for  analysis.  The  torsional  modes  for  N  =  0  were  not 
retained  to  eliminate  a  rotational  input  at  the  base  of  the 
substructure  model.  Modal  contributions  corresponding  to  N 
=  2,3  normal,  meridional  and  circumferential  displacements 
were  retained  for  this  full  model  case.  Figure  2  shows  the 
basic  shapes  of  the  waves  on  the  cylindrical  cross  section. 

C.  INTERNAL  SUBSTRUCTURE  MODELING 

In  ELSHOK  calculations,  a  substructure  model  is  formu¬ 
lated  to  determine  the  modes,  masses  and  natural  frequencies 


of  the  internal  equipment  being  analyzed.  Phase  II  of 
ELSHOK  utilizes  the  SAP  IV  finite  element  code  to  accomplish 
modeling  of  submarine-installed  equipment.  ELSHOK  is 
written  in  program  modules  that  execute  independently  and 
exchange  data  through  the  output  files.  This  implementation 
lends  itself  well  to  the  type  of  parametric  study  conducted 
in  this  thesis  work.  A  variety  of  different  internal  models 
may  be  combined  with  the  same  shell  model  to  gain  insight  on 
the  interaction  of  the  response  between  the  shell  and 
various  substructure  characteristics.  Modeling  the 
substructures  in  a  finite  element  method  more  accurately 
predicts  effects  of  the  substructure  on  the  hull  to  which  it 
is  attached  than  if  the  model  is  represented  as  a 
concentrated  mass.  A  concentrated  mass  model  would  not 
account  for  any  substructure  responses  that  might  become 
dynamically  amplified. 

The  substructure  modeled  for  this  study  is  a  doubly 
overhung  steel  beam  mounted  with  steel  brackets  to  the  port 
side  of  the  shell.  Connection  points  are  at  the  two  large 
rings  at  frames  18  and  27,  placing  the  beam  symmetr ically 
about  the  midships  section.  Only  half  of  the  beam  structure 
requires  modeling  since  symmetry  may  be  employed.  Symmetry 
is  indicated  in  SAP  IV  by  specifying  the  appropriate 
boundary  condition  at  the  point  of  symmetry  on  the  model. 
This  is  similar  to  the  procedure  in  B0S0R4.  The  dimensions 


and  details  of  the  model  are  shown  in  Figure  7.  Figure  7 
also  shows  the  orientation  of  the  substructure  within  the 
shell  and  to  the  shockwave.  The  doubly  overhung  beam  is 
represented  using  the  three-dimensional  beam  element  model 
of  SAP  IV  depicted  in  Figure  8  which  shows  half  the 
structure  with  the  first  mode  for  D  =  4.2  in.  superimposed 
on  it.  In  SAP  IV  nodes  are  used  to  divide  the  model  into 
separate  elements  and  serve  the  purpose  of  call  out  points 
for  velocity-time  histories  in  Phases  III  and  IV  of  ELSHOK. 
Modeling  one  half  of  the  doubly  overhung  beam  was 
accomplished  using  21  nodes,  4  on  the  bracket  and  17  on  the 
beam  half  length.  The  brackets  are  represented  as  purely 
extensional  members  (a  spring  with  mass)  and  the  intro¬ 
duction  of  a  bending  moment  at  the  connection  between  the 
beam  and  bracket  is  eliminated.  This  is  accomplished  by 
specifying  an  end  release  node  at  node  12  that  represents 
the  connection  point  on  the  model.  The  interest  in  the 
response  of  the  substructure  is  limited  to  the  z-direction 
and  restrictions  on  the  degrees  of  freedom  of  the  model  were 
appropriately  specified.  This  reduces  the  calculations 
required  and  produces  responses  in  the  direction  of 
interest . 

In  order  to  investigate  the  response  interaction  between 
the  shell  and  substructures  of  various  weights  and  charac¬ 
teristics  the  depth  D  of  the  beam  is  changed  for  each  case. 


It  should  be  noted  that  the  width  of  the  bracket  varies  with 
the  depth  of  the  beam  by  1.5  x  D.  Developing  the  various 
sized  models  was  accomplished  by  editing  the  SAP  IV  input 
code  to  reflect  the  new  cross-sectional  areas  and  moments  of 
inertia  for  the  beam  and  bracket  elements.  Other  properties 
remained  unchanged. 

D.  SHELL/SUBSTRUCTURE  RESPONSE  MODELING 

The  structural  problem  of  ELSHOK  is  separated  into  two 
parts  to  avoid  calculations  involving  the  modes,  natural 
frequencies  and  stiffness  of  a  combined  system  of  the  shell 
and  substructure.  Even  though  separate  modal  analysis  is 
carried  out  on  each  model  in  different  phases  of  ELSHOK  the 
responses  determined  for  the  shell  and  substructure  are 
arrived  at  by  coupling  the  free-free  modes  of  the  shell  with 
the  fixed-base  natural  modes  of  the  substructure.  In  effect 
the  shell  and  substructure  can  be  approximately  represented 
by  a  two-degrees  of  freedom  mass-spring  model  as  in  Figure 
9.  Here  the  mass  of  the  shell  is  M2  and  the  stiffness 
characteristics  are  K2.  The  substructure  model  is 
represented  by  Ml  and  Kl.  It  has  been  shown  for  some 
situations  using  two  single-degree  of  freedom  models  can 
produce  nearly  the  same  responses  as  a  two-degrees  of 
freedom  model  when  the  mass  ratio  is  very  small.  This  in 
general  does  not  hold  for  systems  of  large  mass  ratios  and 
is  apt  to  produce  marked  differences  in  response  between  the 


two  methods  if  the  same  stiffness  characteristics  are  used 
in  both  [Ref.  10].  ELSHOK  results  can  be  viewed  with  more 
confidence  in  predicting  component  responses  because  it 
utilizes  the  multi-degrees  of  freedom  type  model.  It  is 
preferred  in  this  study  presented  due  to  the  large  mass 
ratios  that  are  involved. 


Figure  6.  Ring-Stiffened  Cylindrical  Shell  Model  Used  in  Analysis 


Substructure  Model  Used  in  Analysis 
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IV.  ANALYSIS 


A.  OBJECTIVES  OF  ANALYSIS 

The  objectives  of  the  analysis  procedure  are  to 
investigate  effects  on  component  responses  brought  about  by 
varying  the  mass  and  stiffness  of  the  internal  substructure. 
More  specifically,  it  is  divided  into  two  major  and  one 
minor  task.  First  of  these  is  to  determine  how  heavy  the 
internal  structure  must  be  before  it  significantly  affects 
the  motion  of  the  shell.  Secondly,  the  dynamic  amplifica¬ 
tion  of  the  responses  of  the  substructure  and  its  effects  on 
the  shell  are  investigated.  Lastly,  it  is  of  interest  to 
know  the  magnitude  of  the  forces  generated  in  the  support 
brackets.  An  integrating  scheme  is  used  to  determine 
deflections  and  thereby  estimate  the  forces. 

The  shell  and  substructure  are  subjected  to  side-on 
shock  loading  from  an  UNDEX.  The  UNDEX  is  represented  in 
two  different  manners.  One  is  an  incident  pressure- t ime 
history  simulating  a  tapered  charge  or  nuclear  detonation 
and  the  other  is  an  empirical  exponentially  decaying 
incident  pressure  simulating  a  spherical  conventional  charge 
detonation.  The  standoff  or  range  from  the  origin  of  the 


UNDEX  to  the  closest 
both  the  tapered  and 


point  of  the  target  is  the  same  for 
conventional  cases.  The  weight  of  the 


conventional  charge  is  selected  to  deliver  the  same  impulse 
as  that  of  the  tapered  charge.  First,  the  transient 
responses  of  the  empty  shell  are  determined  for  a  basis  of 


comparison.  Next,  responses  for  a  variety  of  internal 
structure  sizes  are  calculated  for  both  conventional  and 
tapered  charge  cases.  The  results  are  produced  in  the  form 
of  velocity-time  histories  and  are  obtained  for  points 
specified  on  the  shell  and  substructure. 

B.  ANALYSIS  PROCEDURE 

The  shell  natural  frequencies  and  modes  for  wave  numbers 
N  =  0,  1,  2,  and  3  are  found  using  the  BOSOR4  code.  The 
virtual  mass  array  computed  in  ACESNID  is  done  including 
each  wave  N  specified  in  BOSOR4  for  compatibility.  In  order 
to  simplify  the  analysis  and  since  the  brackets  are  designed 
to  resist  movement  out  of  the  horizontal  plane,  the 
torsional  and  rolling  modes  are  eliminated  from  considera¬ 
tion  to  prevent  a  rotational  input  to  the  base  of  the 
substructure.  This  is  accomplished  by  setting  the  variable 
NTORSN  =  0  and  omitting  any  purely  rolling  or  torsional 
modes  for  N  =  0  from  input  in  the  PIFLASH  code.  All  thirty 
modes  for  the  other  values  of  N  are  used.  This  shell  model 
is  used  throughout  the  study  in  combination  with  each 
substructure.  The  B0S0R4  modes  deleted  and  retained  are 
indicated  in  Tables  I  and  II. 


SAP  IV  was  used  to  determine  the  fixed-base  natural 
modes  and  frequencies  as  D  was  varied  on  the  substructure 


model.  The  investigation  begins  with  a  beam  depth  D  of  one 
inch  and  increases  to  twenty  inches.  It  should  be  noted 
that  twenty  inches  is  nearly  physically  impossible  to  fit 
within  the  confines  of  the  shell  and  was  computed  to 
investigate  the  effects  of  a  large  mass  on  the  shell.  Each 
SAP  IV  run  produced  ten  fixed-base  natural  frequencies 
corresponding  to  different  mode  shapes  of  the  doubly  over¬ 
hung  beam.  Table  III  contains  frequencies  computed  for  the 
cases  D  =  1.0,  4.2  and  9.0  in.  Not  all  the  modes  determined 
in  each  case  were  used  in  the  follow-on  calculations  because 
some  of  the  high  frequencies  would  not  allow  convergence  of 
the  integrating  scheme  used  without  changing  the  time  step. 
All  modes  having  frequencies  of  3979  Hz  and  greater  were  not 
used  in  the  analysis.  These  higher  modes  do  not  contribute 
significantly  enough  to  warrant  varying  the  time  step  to 
remain  within  the  convergence  limits. 

Results  from  each  SAP  IV  run  require  processing  by  the 


PICRUST  code  for  preparation  to  merge  the  shell  and 
substructure  data  files  with  the  fluid  data  files.  The 
configuration  of  each  substructure  installation  along  with 
the  modes  to  be  used  in  Phase  III  calculations  are 
specified.  The  substructure  z-axis  is  related  to  the  shell 
coordinate  system  by  the  angle  shown  in  Figure  5.  The 
meridional  locations  of  the  attachment  points  are  also 
specified  here.  Each  case  investigated  placed  the 


% 
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substructure  attached  to  the  port  side  or  an  angle  of  180 
from  the  shell  Z-axis  orientation  (exactly  one  half  the 
distance  between  the  top  and  bottom) .  The  brackets  connect 
to  the  rings  at  frames  18  and  27  placing  the  substructure 
evenly  about  the  midships  section  of  the  shell  or  node  2/27. 
Each  case  D  is  subjected  to  two  types  of  UNDEX  representa¬ 
tions  having  origins  at  70  ft.  (840  in.)  from  the  shell 
abeam  the  starboard  side  at  midships.  The  origin  also  is 
located  in  the  shell's  Z-X  plane.  After  PICRUST  prepares 
the  shell/substructure-fluid  data  the  time  integration  of 
the  incident  pressures  experienced  by  the  shell  surface  is 
performed  to  determine  velocities. 

The  US LOB  code  models  the  type  of  UNDEX  by  differently 
representing  the  incident  pressures.  The  tapered  charge 
pressure-time  history  is  represented  by  inputting  discrete 
pressure-time  points  that  describe  it.  The  conventional 
spherical  charge  is  represented  by  equation  (1)  and 
specifying  the  appropriate  constants  for  the  type  of 
explosive  used.  The  weight  of  the  charge  specifies  the  size 
of  the  explosion  and  its  resultant  impulse  and  incident 
pressure  on  the  shell.  In  determining  a  charge  weight  that 
would  deliver  the  same  impulse  as  the  tapered  charge 
pressure- time  history  equations  from  Cole  [Ref.  4]  and 
constants  describing  HBX-1  were  used.  Describing  the 
impulse  of  unit  area  of  the  shock  wave  front  up  to  a  time  t 
after  its  arrival  at  the  shell  wet  surface  by: 


I  (t) 


(3) 


^  P(t)dt 

0 

where:  P(t)  =  PQ  exp(t/0Q) 


(4) 


Then  integrating  equation  (2)  over  the  time  after  arrival  of 
the  shock  wave  results  in: 


U t)  =  PO0O[1  -  exp(-t/60)]  (5) 

Substituting  equation  (2)  for  8q  and  PQ  where, 

i  /■»  ^2  • 

PQ  =  K1(Wi/J/R)  (6) 

into  equation  (5)  then  setting  it  equal  to  the  area  under 
the  curve  in  Figure  3  the  charge  weight  arrived  at  through 
iteration  is  352  lbm  of  HBX-1.  The  constants  for  HBX-1 
are  determined  from  data  recorded  in  feet  and  msec.  This 
required  and  K2  to  be  adjusted  to  represent  the  empirical 
relationships  in  inches  and  seconds  in  order  to  be  con¬ 
sistent  with  the  units  used  in  this  application  of  ELSHOK. 

The  standoff  and  lateral  placement  of  the  origin  of  the 
UNDEX  is  specified  in  USLOB  along  with  the  duration  of 
integration  and  the  time  step  to  be  used.  A  time  duration 
of  10  msec  and  a  step  of  3x10“^  sec  was  used  to  capture  the 
peak  response  in  each  case.  Also  no  effects  due  to  surface 


cutoff  or  reflection  of  shock  waves  were  considered.  This 
situation  represents  a  deep  submerged  test  simulation  where 
bulk  cavition  due  to  the  surface  doesn't  occur  and  there  is 
no  bottom  bounce  or  surface  cutoff  effects. 

C.  PLOTTING  PROCEDURE 

Plotting  of  velocity  profiles  is  accomplished  by  PUSLOB 
in  the  last  phase  of  ELSHOK.  Each  profile  is  a  record  of 
the  velocity  history  for  a  single  node  specified  on  the 
shell  or  substructure.  The  plot  is  limited  to  two  angular 
locations  of  the  shell  nodes  and  velocities  are  expressed  in 
in/sec  due  to  the  units  used  in  this  study.  The  version  of 
ELSHOK  available  at  the  Naval  Postgraduate  School  is 
installed  on  the  VAX/VMS  system  and  plots  were  produced  on 
Tektronix  devices.  In  order  to  aid  in  comparing  the 
velocities  of  different  nodes  the  data  files  used  to  create 
the  PUSLOB  plots  were  processed  using  program  PLOTCONV. 
Program  PLOTCONV  rearranges  the  data  files  and  converts  the 
units  of  velocity  from  in/sec  to  ft/sec.  This  facilitates 
the  use  of  plotting  routines  EASYPLOT  or  DISSPLA  available 
on  the  school's  IBM  system  and  expresses  the  plots  in  the 
customary  units  of  the  shock  community:  velocity  in  ft/sec, 
time  in  msec.  Program  FISH1  is  a  program  developed  for  use 
with  the  DISSPLA  plotting  routine.  Copies  of  programs 
PLOTCONV  and  FISH1  are  found  in  Appendix  B. 


D.  FORCES  IN  THE  SUBSTRUCTURE  SUPPORTS 


Since  ELSHOK  produces  results  in  velocity-time 
histories,  force  calculations  require  some  post  processing. 
Forces  generated  in  the  support  brackets  would  be  mainly  due 
to  extension  and  compression  because  the  bracket  design 
resists  motion  out  of  the  horizontal  plane.  Due  to  this  and 
to  simplify  the  analysis,  inertia  of  the  brackets  is 
neglected  in  the  force  calculations,  though  their  mass  is 
included  in  the  shock  response  calculations  of  ELSHOK. 
Program  FORCES  found  in  Appendix  C  was  used  to  convert  the 
velocity  histories  to  deflections  using  a  Simpson's  one- 
third  integrating  scheme.  Uniaxial  forces  estimated  are 
based  on  F  =  eEA  where  e  is  the  estimated  deflection  of  the 


bracket  divided  by  the  length  between  the  connection  points 


of  bracket  at  the  ring  and  beam  junction.  E  is  the  modulus 
of  elasticity  of  steel  and  A  is  the  cross  section  area  of 


the  bracket,  1.5  x  D. 
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TABLE  III 


Mode 


=  1.0  in . 


Frequency  (Hz) 


D  =  4.2  in . 


0  =  9.0  in 


25.46 

106.9 

228.7 

81  .63 

330.5 

616 .3 

337.4 

1354.0 

1728  .0 

463.0 

1421.0 

3032.0 

966.4 

2700.0 

4982.0 

1206.0 

4167.0 

6552  .0 

1802.0 

5475.0 

8951  .0 

2151  .0 

6728  .0 

9957.0 

2680.0 

8380.0 

12430.0 

3281  .0 

9939  .0 

17810.0 

V.  RESULTS 


Numerous  cases  were  explored  in  this  parametric  study  in 
which  the  beam  depth  D  varied  over  the  range  of  one  to 
twenty  inches.  ELSHOK  generated  a  large  amount  of 
information  about  each  case.  The  most  signiiicant  results 
have  been  summarized  in  tabular  and  plotted  forms  which 
contain  information  on  cases  subjected  to  the  tapered  charge 
and  the  conventional  charge  UNDEX  simulations.  SAP  IV 
predicted  frequencies  and  modes  for  the  tapered  charge  case 
with  beam  depth  D  =  4.0  in.  that  compared  very  closely  with 
the  ONR  model  previously  tested.  This  provides  a  certain 
amount  of  confidence  in  the  SAP  IV  model  used  throughout 
this  investigation. 

A.  SHELL  AND  SUBSTRUCTURE  INTERACTION 

Even  though  the  variation  of  substructure  sizes  included 
the  cases  resulting  near  the  600  Hz  secondary  frequency  that 
was  expected  to  induce  dynamic  amplification  in  the  original 
ONR  model  there  is  no  such  response  over  the  entire  range  of 
cases  investigated.  The  case  for  D  =  4.2  in.  conventional 
charge  resulted  in  a  peak  of  velocities  at  the  tip  and 
midsection  of  the  beam  as  shown  in  Figure  10  and  Figure  11. 
However,  this  is  not  a  resonant  condition.  This  slight 
amplification  can  be  traced  back  to  shell  modes  1,  2,  3  for 


N  =  2  and  mode  1  for  N  =  3.  Tables  IV  and  V  summarize  the 
range  of  results  for  both  charge  types  and,  when  viewed  with 
Figures  12  and  13,  there  is  some  indication  of  transfer  of 
energy  between  the  substructures  and  the  hull,  but  not  much. 
There  seems  to  be  no  significant  interaction  between  the 
beam  and  shell.  The  basic  shape  of  the  velocity-time 
histories  at  the  beam  tip,  midsection  and  junction  remain 
nearly  the  same  throughout  the  range  of  D  studied.  The  only 
effects  seem  to  be  the  decrease  of  beam  and  bracket  maximum 
velocities  as  the  weight  of  the  beam  is  increased.  The 
complete  velocity  profiles  for  the  beam  bracket  and  connec¬ 
tion  points  for  the  D  =  4.2  in.  cases  are  shown  in  Figures 
14  through  17.  This  size  was  selected  for  comparison 
because  it  did  demonstrate  slightly  higher  overall  velo¬ 
cities  on  the  beam  and  also  falls  nearly  into  mid-range  of 
the  cases  if  D  =  15.0  and  20.0  in.  are  regarded  as  unreason¬ 
able  models  due  to  their  size.  These  two  cases  were  looked 
at  simply  to  check  the  effect  of  the  large  mass  they  would 
represent  on  the  hull.  Noting  Figure  13  that  shows  the 
maximum  velocities  of  the  bracket  connection  points,  the 
rings  at  frames  18  and  27,  it  is  seen  that  the  variation 
in  peak  responses  is  only  slight  over  the  range  of  D  = 
4.0-10.0  in.  for  the  tapered  charge  case.  This  corresponds 
to  substructure  to  shell  mass  ratios  of  .07644- .  19875  .  The 
empty  shell  velocity  at  this  point  is  10.48  ft/sec  and  its 


highest  maximum  velocity  occurs  when  D  =  7.0  in.  which  is 
11.88  ft/sec  corresponding  to  a  mass  ratio  of  .13638.  The 
conventional  charge  empty  shell  velocity  at  this  point  was 
32.14  ft/sec  and  this  was  the  maximum  velocity  recorded  on 
frames  18  and  27  of  the  conventional  cases. 

The  velocity-time  histories  of  the  beam  and  bracket  for 
the  cases  D  =  2.0-10.0  in.  for  both  charge  types  are  very 
similar  to  the  D  =  4.2  in.  cases  shown  in  Figures  14  through 
16.  At  D  =  1.0  in.  the  response  was  very  oscillatory  and 
produced  the  largest  maximum  velocity  for  the  conventional 
charge  at  the  beam  tip.  This  could  be  due  to  the  large 
impulse  combined  with  the  smallest  mass  ratio  of  .01874. 

The  other  cases  simply  show  overall  decreases  in  velocities 
but  retained  the  basic  shape  of  the  velocity-time  histories. 
Comparing  the  empty  shell  to  the  shell  with  substructure 
velocity-time  histories  at  the  midsection  of  the  shell  and 
frame  18  demonstrates  what  little  contribution  the 
substructures'  responses  made  in  the  shell  response.  In 
Figures  17  through  20  the  shape  of  the  velocity-time 
history  profiles  for  these  points  can  be  seen  to  be  mostly 
unchanged  by  the  presence  of  the  substructure.  This  was 
true  for  all  cases,  conventional  and  tapered,  and  increasing 
the  weight  of  the  substructure  served  mostly  to  reduce  the 
velocity  at  the  connection  point. 

Comparing  the  frequencies  of  the  shell  model  from  Table 
II  with  the  first  five  natural  frequencies  of  the  substruc- 


ture  models  there  are  very  few  matched  frequencies  at  and 
below  600  Hz  which  are  the  range  expected  to  excite  the 
first  two  modes  of  the  substructure.  The  frequencies  for 
the  cases  listed  in  Table  III  show  how  quickly  they  increase 
to  magnitudes  where  their  contributions  to  response  is 
minimal.  Any  modes  of  the  substructure  being  excited  by  the 
base  input  at  the  connection  points  of  the  brackets  and 
shell  did  not  result  in  a  dynamically  amplified  response. 
Each  increment  of  substructure  size  produces  a  set  of  fixed 
base  frequencies  of  larger  magnitude  than  the  previous  set. 
If  the  natural  frequency  is  equal  to: 


where  K  =  stiffness  and  =  effective  mass  then  the 
stiffness  must  be  increasing  more  rapidly  than  the  mass  of 
each  structure  to  produce  the  higher  frequencies. 

B.  OVERALL  SHELL  RESPONSES 

The  shell  response  overall  seems  to  exhibit  some  of  the 
expected  modes  for  the  circumferential  waves  used  in  the 
calculations  with  the  exception  of  the  purely  torsional 
modes  that  were  deleted  early  on.  First  looking  at  the 
velocity-time  histories  in  Figures  21  through  24  the  extreme 
difference  in  velocities  produced  by  the  two  types  of 
charges  is  apparent.  As  expected,  the  conventional  charge 


produces  larger  velocities  in  each  case.  These  velocities 
are  in  the  athwart  ship  direction,  positive  sense  to  port, 
at  call  out  points  on  the  midships  station  and  the  three 
frames  where  the  discrete  rings  are  located.  The  velocity¬ 
time  histories  due  to  the  conventional  charge  illustrate  how 
the  velocities  propagate  from  the  midships  station  to  the 
end  plates.  Looking  at  the  peak  velocities  the  shell  can  be 
seen  to  be  bending  about  the  midships  station  while  the 
entire  shell  translates  in  the  direction  of  the  shock  wave. 
These  modes  can  be  seen  in  the  velocity  profiles  for 
conventional  and  tapered  charge  cases  where  the  midships 
peak  velocity  is  less  than  those  for  the  two  rings  located 
at  frames  9  and  18.  Looking  at  Figures  22  and  25  shows  that 

the  velocities  at  the  midships  station  for  the  conventional 
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case  on  the  port  and  starboard  sides  are  of  opposite  signs 
at  about  time  equal  to  1. 0-2.0  msec.  This  illustrates  a 
higher  order  circumferential  wave  of  N  =  2.  Figures  25  and 
26  show  the  velocity  profiles  at  the  top  of  the  shell  in  the 
athwart  ship  direction.  In  each  case  run  the  top  and  bottom 
velocity  profiles  were  exactly  matched  indicating  the 
symmetric  response  to  the  shock  loading  in  the  translational 
and  whipping  modes.  Figures  27  and  28  illustrate  the  slight 
lag  in  the  velocities  on  the  port  side.  As  a  rule  the 


conventional  starboard  side  velocities  at  ring  frames  were 
always  higher  than  t-he  po.'t.  This  is  due  to  the  short 


m 


duration  of  shock  front  of  the  conventional  charge. 
Alternately,  the  tapered  charge  caused  slightly  higher  peak 
velocities  on  the  port  side  but  also  lagged  the  peak 
responses  on  the  stbd  side.  This  shows  how  the  tapered 
charge  allows  the  port  side  to  exhibit  more  response  to  the 
exciting  forces. 

The  end  plate  velocities  are  shown  in  Figures  29  and  30. 
The  end  plate  in  Figure  29  can  be  seen  to  bulge  in  and  out 
at  the  center  for  the  conventional  case  and  in  Figure  30  the 
tapered  case  shows  the  same  response.  The  ring  at  frame 
zero  tends  to  hold  the  outer  circumference  of  the  end  plate 
in  place  as  shown  by  the  lower  velocities  experienced  at  the 
outer  edge.  These  velocity-time  histories  are  indicative  of 
breathing  modes  of  the  shell. 

C.  FORCES  IN  THE  SUPPORTS 

The  axial  forces  were  estimated  for  the  D  =  4.2  in. 
model  subjected  to  both  the  conventional  and  tapered 
charges.  In  both  cases  the  forces  fluctuate  between 
compressive  and  tensile  as  indicated  by  the  changes  in  sign 
in  the  results  in  Appendix  C.  The  maximum  force  estimated 
occurred  in  the  conventional  case.  It  was  a  compressive 
force  of  about  103  kips  or  a  stress  of  about  16  psi.  The 
maximum  tensile  force  experienced  by  the  support  bracket  for 
this  case  was  estimated  at  approximately  75  kips.  The 
tapered  charge  case  produced  much  lower  force  estimates. 


The  maximum  was  tensile  and  occurred  at  about  1.7  msec 
whereas  the  conventional  maximum  force  occurred  at  about 
.3  msec.  The  maximum  estimated  tensile  and  compressive 
forces  for  the  tapered  charge  case  are  25  kips  and  23  kips 
respectively . 


D(in)  Fl(Hz)  F2 ( H  z ) 


MAXIMUM  VELOCITIES  (FPS) 


M1/M2 


TIP( 04 ) 


1 

25.46 

81  .63 

27.92 

2 

50.92 

162.2 

24.58 

3 

76.38 

240.6 

22.10 

4 

101  .8 

315.9 

21.18 

4.2 

106.9 

330.5 

20.71 

5 

127.3 

387.1 

18.98 

6 

152.8 

453.4 

17.49 

7 

1 7  8 . 2 

514.0 

15.69 

8 

203.7 

568.7 

13.72 

9 

228.7 

616 .3 

12.05 

10 

254.6 

659.6 

12.04 

15 

381  .9 

798.5 

11.16 

20 

509.2 

863,3 

10.79 

(Shell  mass  M2  =  3420.0  lhm, 
Fnpty  shell  V(2/21)  =  10.48 


MID( 2 1 )  JCN (12)  2/21(180°) 


19.76 

10.96 

10.32 

.01874 

21  .79 

11.39 

10.19 

.0377 

19.14 

11  .02 

10.11 

.05696 

18.85 

11.94 

11.37 

.07644 

18.52 

11  .97 

1 1  .53 

.0765 

17.21 

11.98 

11  .76 

.09602 

15.87 

11  .89 

11  .72 

.11615 

14.33 

12.12 

1 1  .88 

.13638 

12.70 

11  .83 

11  .69 

.15685 

11  .76 

11.43 

11  .24 

.1775 

1 1  .75 

11.11 

10 .9« 

.19875 

11.17 

10.74 

10.81 

.3071 

10.75 

10.68 

10.37 

.4219 

FPS) 

£ 


TABLE  V 

SUMMARY  OF  RESULTS — CONVENTIONAL  CHARGES 


D(in)  FI ( Hz)  F2 ( Hz )  MAXIMUM  VELOCITIES  (FPS) 

TIP ( 04  )  MID( 21 )  J CN (12)  2/21(180°) 


1 

25.46 

81  .63 

77.21 

43.93 

33  .97 

31  .75 

.01874 

2 

50.92 

162.2 

55.66 

41  .32 

31  .97 

29.81 

.0377 

3 

76.38 

240.6 

51  .59 

39.90 

33.44 

27.85 

.05696 

4 

101  .8 

315.9 

58.51 

48.48 

30.77 

28.02 

.07644 

4.2 

106.9 

330.5 

59.16 

48.91 

29.87 

27.83 

.0765 

5 

127.3 

387.1 

58.04 

48.76 

26.56 

26.43 

.09602 

6 

152.8 

453.4 

54.42 

47.06 

24.18 

23.45 

.11615 

r 

7 

178.2 

514.0 

49.22 

43.05 

24.48 

21.78 

.13638 

8 

203.7 

568.7 

43.33 

38.35 

25.06 

23  .02 

.15685 

9 

228.7 

616.3 

37.78 

33.81 

34.55 

23.33 

.1775 

10 

254.6 

659.6 

32.35 

29.65 

23.48 

22.90 

.19875 

15 

381  .9 

798.5 

20.91 

20.34 

19.65 

19.69 

.3071 

20 

509.2 

863.3 

17.44 

17.47 

17.12 

16.95 

.4219 

(Shell  mass  M2  =  3420.0  lbm. 

Empty  shell  V ( 2/2 1 )  =  34.14  FPS) 


MAX  VELOCITIES  VS  BEAM  DEPTH 


Figure  10.  Maximum  Velocities  at  Tip  of  Beam 
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VELOCITY  (FT/SEC) 


VELOCITY  PROFILE- 


MID (21 )  V-flTH 


Figure  15  . 


Velocity-Time  Histories  at  Beam  Midsection  for 
D  =  4 . 2  In. 
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VELOCITY  (FT/SEC) 


VELOCITY  PROFILE 


JCN( 12)  V-flTH 


Figure  16.  Velocity-Time  Histories  at  Beam/Bracket  Junction 
for  D  =  4.2  In. 
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Figure  18.  Velocity-Time  Histories  at  Frame  18  for  Emotv 
Shell 
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VELOCITY  l FT/SEC) 


VELOCITY  (FT/SEC) 


VELOCITY  PROFILE 

NODE  2/27  V-flTH  PORT 


TIME  (MSEC) 


Figure  21.  Velocity-Time  Histories  at  Midships  with 
Substructure  of  D  =  4.2  In. 
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VELOCITY  (FT/SEC) 


VELOCITY  PROFILE 


NODE  2/27  V-flTH  PORT 


re  22.  Velocity-Time  Histories  at  Midshins  for  Emoty 
Shell 
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VELOCITY  PROFILE 


NODE  2/12  V-flTH  PORT 


Figure  23.  Velocity-Time  Histories  at  Frame  9  for  F.moty 
Shell 
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END  PLATE  CONV  V-FWD 


tre  29.  Velocity-Time  Histories  at  End  Plate  for 
Conventional  Charge 
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END  PLATE  TAP  V-FWD 


Figure  30.  Velocity-Tine  Histories  at  End  Plate  for  Tapered 
Charge 


VI.  CONCLUSIONS 


The  results  from  the  cases  studied  tend  to  show  that  a 
model  can  be  designed  to  minimize  dynamic  amplification. 

In  no  case  did  a  truly  resonant  situation  occur.  As  the 
substructure  model  increased  in  size  and  mass  so  did  the 
stiffness  and  the  effect  was  as  if  a  lumped  mass  was 
attached  to  the  shell.  Even  though  substructure  sizes 
having  frequencies  that  matched  those  of  the  hull  were 
tested,  the  corresponding  mode  shapes  being  excited  did  not 
result  in  dynamic  amplification  of  the  beam  responses.  The 
beam  model  having  a  secondary  frequency  of  330.5  Hz  became 
slightly  amplified  but  only  for  the  conventional  explosive 
analysis . 

Though  the  models  used  in  this  study  are  based  on  real 
configurations  it  is  important  to  keep  in  mind  this  was  a 
numerical  study  on  a  simplified  model  of  an  already  simple 
structure.  Even  though  the  substructure  design  used  did  not 
become  highly  excited  when  combined  with  the  shell  in  this 
study,  the  beam’s  behavior  might  change  if  combined  with  a 
shell  that  has  a  frequency  range  and  mode  shapes  that  are 
different  from  those  of  the  shell  used  in  this  study.  The 
shell  model  was  a  relatively  small  structure  that  resulted 
in  a  majority  of  frequencies  too  high  to  excite  the 
substructure  significantly. 
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ELSHOK  has  the  capability  of  representing  much  more 
complex  models  than  those  used  here.  It  has  the  capability 
of  simulating  a  full  sized  submarine  hull  with  internally 
attached  equipment.  The  ELSHOK  code  has  been  shown  to  give 
fairly  accurate  results,  especially  in  the  early  portion  of 
shock  problems,  through  comparisons  of  its  predictions  to 
physical  test  results.  However,  like  any  numerical  device 
it  needs  such  comparisons  to  verify  its  results.  On  the 
other  hand,  conducting  tests  using  ELSHOK  has  the  distinct 
advantage  of  eliminating  risk  of  damage  to  valuable  assets. 
Once  a  model  has  been  verified  to  satisfactorily  predict  the 
prototype  responses,  a  variety  of  internal  equipment  models 
could  be  applied  to  perform  a  variety  of  studies  at 
fractions  of  the  cost  of  building  real  models.  Again  the 
drawbacks  are  the  pitfalls  of  any  numerical  endeavor, 
namely,  the  output  is  only  as  good  as  the  input.  Simplified 
models  may  produce  simplified  results  and  poor  assumptions 
must  be  guarded  against. 

It  is  important  to  realize  that  the  results  of  this 
study  were  arrived  at  through  the  implementation  of  a 
potentially  powerful  tool  for  future  research  at  the  Naval 
Postgraduate  School,  namely  ELSHOK.  Further  studies  in  the 
area  of  UNDEX  using  the  code  should  be  carried  out  to  prove 
its  validity  or  develop  it  into  more  useful  applications  of 
design  and  research.  Suggestions  on  improvements  to  develop 


the  line  of  research  in  this  thesis  from  a  numerical 
sensitivity  study  into  a  realistic  problem  are  to  obtain  or 
construct  a  representative  full-size  submarine  model  or  SSTV 
for  B0S0R4  and  model  all  internal  equipments  to  be  tested  in 
appropriate  detail  while  applying  fewer  restrictions  to  the 
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APPENDIX  A 

INPUT  FILES  FOR  ELSHOK 

1.  BOSOR4  INPUT  DATA  FOR  ALL  CASES 

The  following  data  file  was  used  to  generate  the 
cylindrical  shell  model.  Four  separate  runs  of  B0S0R4  were 
made  corresponding  to  N  =  0,1, 2, 3  by  setting  NOB,  NMINB  and 
NMAXB  equal  to  each  case  N.  All  other  inputs  remain  the 


same  for  each  case 


RING- STIFFENED 


16.  812^0 

J 

2. 500000 

N  0 
0 

8 

2 

0. 1080000E+08 
0. 3200000 
0. 2535000E-03 

8 


16.81250 

0.  1875000 
N 


3. 500000 
53. 12500 
103. 7500 


CLINDER  WITH  FLAT  ENDPLATES-- SYMMETRY  USED--N=Q 
INDIC  3  ANALYSIS  TYPE  INDICATOR  „  „  v 

NPRT  3  OUTPUT  OPTIONS  (^MINIMUM.  2=MEDIUM.  33MAXIMUM)  v 
ISTRES3  OUTPUT  CONTROL  (^RESULTANTS,  1=SIGMA.  2=EPSIL0N) 
NSEG  *  NUMBER  OF  SHELL  SEGMENTS  (LESS  THAN  25) 


SEGMENT  NUMBER  1  1 

NMESH  3  NUMBER  OF  NODE  P' 
NTYPEH3  CONTROL  INTEGER 
NSHAPE3  INDICATOR  (1.2  0: 
R1  3  RADIUS  AT  BtGIN 


^S0R  ^F^NODAL^oi^TT  ^ 


OF  DISCRE' 
)  FOR  TYPE 


F.  TO  REFERENCE  SURF. 

S).  ETC.  FOR  THIS  SEGMENT? 
E  RINGS  IN  THIS  SEGMENT 
LB/IN**3_)IN  THIS  SEG. 

OR  TYPE  OF  LINE  LOADS 
OF  SURFACE  LOADING 
FOR  WALL  CONSTRUCTION 


0. 3000000E-08 
7.  125000 
7. 535500 
2. 375000 

-1.  781250 
0 

. 7501E*08 
0. 7330000E-03 
0. 3000000E*08 
8.250000 
11.  69800 
2.  750000 

-2. 062500 

0. 9249000K*08 
0.  7330Q00F.-03 
0.  3000000E.+08 
8. 250000 


PACING 


M*  RADIUS  AT  END  OF  SEGMENT 

3  AXIAL  COORDINATE  AT  END  OF  SEGMENT 
IMP  *  INDICATOR  FOR  IMPERFECTION  (0=NONE,  1=S0ME) 

NTYPEZ3  CONTROL  (1  OR  3 )  FOR  REFERENCE  SURFACE  LOCATION 
ZVAL  3  DISTANCE  FROM  LEFTMOST  SURF.  TO  REFERENCE  SURF. 

DO  YOU  WANT  TO  PRINT  OUT  R(S).  RY?s).  ETC.  FOR  THIS  SEGMENT? 
NRINGS3  NUMBER  (MAX=20J  OF  DISCRETE  RINGS' IN  THIS  SEGMENT 
K= ELASTIC  FOUNDATION  MODULUS  (E.G.  LB /IN** 3 UN  THIS  SEG. 
LINTYP3  INDICATOR  (0.  1,  2  OR  3)  FOR  TYPE  OF  LINE  LOADS 
NLTYPE= CONTROL  (0,1,2.35  FOR  TYPE  OF  SURFACE  LOADING 
NWALL= INDEX  (1.2!  4.  3;  6.  7.  8J  FOR  WALL  CONSTRUCTION 
E  3  YOUNG’S  MODULUS  FOR  SKIN 

U  3  POISSON'S  RATIO  FOR  SKIN  _  .  ___  _  „ 

SM  3 MASS  DENSITY  OF  SKIN  Tz. G.  ALUM. 3. 00025  LB-SEC**2/IN**4) 
ALPHA  3  COEFFICIENT  OF  THERMAL  EXPANSION 

ANRS  3  CONTROL  (0  OR  1)_F0R  ADDITION  OF  SMEARED  STIFFENERS 
SUR  3  CONTROL  FOR  THICKNESS  INPUT  (0  OR  1  OR  -1) 

DO  YOU  WANT  TO  PRINT  OUT  THE  Cft.J)  AT  MERIDIONAL  STATIONS? 

DO  YOU  WANT  TO  PRINT  OUT  DISTRIBUTED  LOADS  ALONG  MERIDIAN? 

SEGMENT  NUMBER  2  2  2,2  2  2  2  %  2 

NMESH  3  NUMBER  OF  NODE  POINTS  (5  3  MIN.;  98  3  MAX. )(_2) 
NTYPEH3  CONTROL  INTEGER  [ l  OR  2  OR  3 J  FOR  NODAL  POINT  SPACING 
NSHAPE3  INDICATOR  (1,2  OR  4]  FOR  GEOMETRY  OF  MERIDIAN 
R1  3  RADIUS  AT  BEGINNING  OF  SEGMENT  (SEE  P.  66) 

Z1  3  AXIAL  COORDINATE  AT  BEGINNING  OF  SEGMENT 

R2  3  RADIUS  AT  END  OF  SEGMENT 

Z2  3  AXIAL  COORDINATE  AT  END  OF  SEGMENT 

IMP  3  INDICATOR  FOR  IMPERFECTION  (0=NONE.  13S0ME5 

NTYPEZ3  CONTROL  ( 1  OR  3)  FOR  REFERENCE  SURFACE  LOCATION 

ZVAL  3  DISTANCE  FROM  LEFTMOST  SURF.  TO  REFERENCE  SURF.. 

DO  YOU  WANT  TO  PRINT  OUT  R(S).  RTS).  ETC.  FOR  THIS  SEGMENT? 
NRINGS3  NUMBER  (MAX320)  OF  DISCRETE  AlNGS  IN  THIS  SEGMENT 
NTYPE  3  CONTROL  FOR  IDENTIFICATION  OF  RING  LOCATION  (2=Z,  3=R 

Z(I)  3  AXIAL  COORDINATE  OF  ITH  RING,  zf  1) 

Z(I)  3  AXIAL  COORDINATE  OF  ITH  RING,  ZC  2) 

zm  3  AXIAL  COORDINATE  OF  ITH  RING,  Zf  55  ... 

NTYPER3  TYPE  ( 0  OR  1  OR  2  OR  4  OR  5 )  OF  DISCRETE  RING  NO. (  1) 

NTYPER3  TYPE  ? 0  OR  1  OR  2  OR  4  OR  5 )  OF  DISCRETE  RING  NO.  f  ZJ 

NTYPER3  TYPE  f 0  OR  1  OR  2  OR  4  OR  5 )  OF  DISCRETE  RING  NO. (  3) 

E  3  young's  MODULUS  OF  RINGf  1)  . 


ANT  TO  PRINT  OUT  R(S).  RT(S),  ETC.  FOR  THIS  SEGME 
NUMBER  (MAX3 20)  OF  DISCRETE  AlNGS  IN  THIS  SEGMENT 
CONTROL  FOR  IDENTIFICATION  OF  RING  LOCATION  (2=Z,  _  ... 
AXIAL  COORDINATE  OF  ITH  RING,  ZC  1) 

AXIAL  COORDINATE  OF  ITH  RING  ZC  2) 

AXIAL  COORDINATE  OF  ITH  RING,  Zf  55  ,  ,  , 

TYPE  (0  OR  1  OR  2  OR  4  OR  5 )  OF  DISCRETE  RING  NO. f  1) 

TYPE  (0  OR  1  OR  2  OR  4  OR  5)  OF  DISCRETE  RING  NO.  f  ZJ 

TYPE  (0  OR  1  OR  2  OR  4  OR  5)  OF  DISCRETE  RING  NO. (  3) 

YOUNG'S  MODULUS  OF  RINGf  1)  „ 

CROSS  SECTION  AREA  OF  RINGf  1)  ,  _  w  ,  „ 

MOMENT  OF  INERTIA  ABOUT  Y-AXIS  (SEE  FIG.  ON  P. 70)(  1) 
MOMENT  OF  INERTIA  ABOUT  X-AXIS(  1) 

PRODUCT  OF  INERTIAL  1]  ,  .  w  .  , 

RADIAL  COMPONENT  OF  RING  ECCENTRICITY, ( SEE  P.  70 )(  1) 

AXIAL  COMPONENT  OF  RING  ECCENTRICITY(  I) 

TORSIONAL  RIGIDITY?  I),  _ 

RING  MATERIAL  DENSITY  LS- G.  ALUMINUM3. 0002535 )(  1) 
YOUNG'S  MODULUS  OF  RINGf  2)  . 

CROSS  SECTION  AREA  OF  RING?  2)  ,  ,, 

MOMENT  OF  INERTIA  ABOUT  Y-AXIS, (SEE  FIG.  ON  P. 70 )  (  2) 
MOMENT  OF  INERTIA  ABOUT  X-AXIS(  2) 

PRODUCT  OF  INERTIAL  2]  ,  .... 

RADIAL  COMPONENT  OF  RING  ECCENTRICITY  (SEE  P.  70 )(  2) 

AXIAL  COMPONENT  OF  RING  ECCENTRICITY(  2) 

TORSIONAL  RIGIDITY?  2) 

RING  MATERIAL  DENSITY  G.  ALUMINUM3. 0002535)(  2) 
YOUNG'S  MODULUS  OF  RINGf  3) 

CROSS  SECTION  AREA  OF  RING?  3) 


OF  DISCRETE  RING  NO 
OF  DISCRETE  RING  NO 
OF  DISCRETE  RING  NO 


:(!} 


70)(  1) 
70)(  1) 


i  V7«  iwi 

2) 

f-AXIS 


ALUMINUM3. 0002535 )(  1) 


11. 69800 
2.  750000 

-2. 062500 

0. 9249000E+08 
0. 7330000E-03 
0 
0 
0 
2 

0. 3000000E+08 
0. 3000000 
0. 7330000E-03 
0 


N 

Y 

0 

3. 3000000E+08 
3. 3000000 
3. 7330000E-03 

Y 

5. 615000 
3.2500000 
1.  250000 
N 
N 
H 
H 

-1 

0 

0 

0 

1 

30 

0 


H 

H 


H 

H 


1 

0 

0 

1 

0 

0 

Y 

Y 


7 

1 

1 

1 

1 

0 

0 

Y 

Y 


IY 

IX 

IXY 

El 

E2 

GJ 

RM 


MOMENT  OF  INERTIA  ABOUT  Y-AXIS, ( S^E  FIG.  ON  P. 70 )( 


3) 


MOMENT  OF  INERTIA  ABOUT  X-AXIS(  3] 

PRODUCT  OF  INERTIAL  3]  , 

RADIAL  COMPONENT  OF  RING  ECCENTRICITY  ( SEE  P.  70 )(  3) 
AXIAL  COMPONENT  OF  RING  ECCENTRICITY(  3) 

TORSIONAL  RIGIDITY?  3) 

RING  MATERIAL  DENSITY  (.E.  G.  ALUMINUM*. 0002535 )(  3) 
-  - -  .S.  -  LB/IN**3)IN  this  seg. 


mis  SEGMENT? 
(Y  OR  N)? 
is  (CONSTANT) 


K= ELASTIC  FOUNDATION  MODULUS  (E.G.  LB/IN**3)IN  THIS  SEG. 
LINTYP*  INDICATOR  (0.  1.  2  OR  3)  FOR  TYPE  OF  LINE  LOADS 
NLTYPE* CONTROL  (0,1, 2, 3)  FOR  TYPE  OF  SURFACE  LOADING 
NWALL* INDEX  fl.  2.  4,  5;  6.  7.  8)  FOR  WALL  CONSTRUCTION 
E  =  young* S  m6du£us  foA  sAin 

U  =  POISSON'S  RATIO  FOR  SKIN 

SM  *MASS  DENSITY  OF  SKIN  (E.G.  ALUM. =. 00025  LB-SEC**2/IN**4) 
ALPHA  =  COEFFICIENT  OF  THERMAL  EXPANSION 

ANRS  =  CONTROL  To  OR  ILFOR  ADDITION  OF  SMEARED  STIFFENERS 
SUR  =  CONTROL  FOR  THICKNESS  INPUT  (0  OR  1  OR  -1) 

ARE  THERE  STRINGERS  (PLEASE  ANSWER  Y  OR  H  ? 

ARE  THERE  RINGS  [ PLEASE  ANSWER  Y  OR  N)? 

K2  = CONTROL  ( 0  OR  1)  FOR  INTERNAL  OR  EXTERNAL  RINGS 
E2  =  RING  MODULUS 

U2  =  RING  POISSON  RATIO 

RGMD=  RING  MASS  DENSITY 
IS  THE  RING  CROSS  SECTION  CONSTANT  IN 
IS  THE  RING  CROSS  SECTION  RECTANGULAR  „ 

D2  =  ARC  LENGTH  BETWEEN  ADJACENT  RING 
T2  =  THICKNESS  OF  RING  (CONSTANT) 

H2  =  HEIGHT  OF  RING  (.CONSTANT) 

DO  YOU  WANT  TO  PRINT  OUT  THE  Cfl.J)  AT  MERIDIONAL  STATIONS? 

DO  YOU  WANT  TO  PRINT  OUT  DISTRIBUTED  LOADS  ALONG  MERIDIAN? 

GLOBAL  DATA _ 

NLAST  =■  PLOT  OPTIONS  f-l=NONE.  0=GEOMETRY.1=U,V.W) 

NOB  =  STARTING  NUMBER  OF  CIRi.  WAVES  (BUtKLINd  ANALYSIS 
NMINB  =  MINIMUM  NUMBER  OF  CIRC.  WAVES  (BUCKLING  ANALYSIS 
NMAXB  =  MAXIMUM  NUMBER  OF  CIRC.  WAVES  (.BUCKLING  ANALYSIS 
INCRB  =  INCREMENT  IN  NUMBER  OF  CIRC.  WAVES  (BUCKLING) 

NVEC  =  NUMBER  OF  EIGENVALUES  FOR  EACH  WAVE  NUMBER 
P  =  PRESSURE  OR  SURFACE  TRACTION  MULTIPLIER 

TEMP  *  TEMPERATURE  RISE  MULTIPLIER 

OMEGA  =  ANGULAR  VEL.  ABOUT  AXIS  OF  REVOLUTION  (RAD/ SEC) 

NUMBER  OF  POLES  (PLACES  WHERE  R=0)  IN  SEGMENT (  1) 

CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO. 

IPOLE  =  NODAL  POINT  NUMBER  OF  POLE,  IPOLEf 

AT  HOW  MANY  STATIONS  IS  THIS  SEGMENT  CONSTRAINED  TO  GROUND? 

IS  THIS  SEGMENT  JOINED  TO  ANY  LOWER- NUMBERED  SEGMENTS? 

NUMBER  OF  POLES  (PLACES  WHERE  R=0)  IN  SEGMENT (  2) 

CONSTRAINT  CONDITIONS  FOR  SEGMENT  NO.  2222 
AT  HOW  MANY  STATIONS  IS  THIS  SEGMENT  CONSTRAINED  TO  GROUND? 
INODE  =  NODAL  POINT  NUMBER  OF  CONSTRAINT  TO  GROUND,  INODE(  1) 
IUSTAR* AXIAL  DISPLACEMENT  CONSTRAINT  ( 0  OR  1  OR  2 ) 

- - -  - - - '  "  2  =  IMPOSED) 

= IMPOSED) 

_  .  .  _  _ •  IMPOSED) 

D1  =  RADIAL  COMPONENT  OF  OFFSET  OF  GROUND  SUPPORT 
D2  =  AXIAL  COMPONENT  OF  OFFSET  OF  GROUND  SUPPORT 
IS  THIS  CONSTRAINT  THE  SAME  FOR  BOTH  PREBUCKLING  AND  BUCKLING? 
IS  THIS  SEGMENT  JOINED  TO  ANY  LOWER- NUMBERED  SEGMENTS? 

AT  HOW  MAY  STATIONS  IS  THIS  SEGMENT  JOINED  TO  PREVIOUS  SEGS. ? 
INODE  =  NODE  IN  CURRENT  SEGMENT  riSEGJ  OF  JUNCTION,  INODEL  l) 
JSEG  =  SEGMENT  NO.  OF  PREVIOUS  SEGMENT  INVOLVED  Ilfr  JUNCTION 
.JNODE  =  NODE  IN  PREVIOUS  SEGMNT  USEG)  OF  JUNCTION 
IUSTAR*  AXIAL  DISPLACEMENT  fO=NOT  SLAVED.  1=SLAVED) 

IVSTAR*  CIRCUMFERENTIAL  DISPLACEMENT  (0=NOT  SLAVED,  1=SLAVED) 
IWSTAR*  RADIAL  DISPLACEMENT  (0=NOT  SLAVED,  l=SLAVEf>) 

I CHI  =  MERIDIONAL  ROTATION  (_0=NOT  SLAVED,  1=SLAVED) 

D1  =  RADIAL  COMPONENT  OF  JUNCTURE  GAP 
D2  =  AXIAL  COMPONENT  OF  JUNCTURE  GAP 

IS  THIS  CONSTRAINT  THE  SAME  FOR  BOTH  PREBUCKLING  AND  BUCKLING? 
DO  YOU  WANT  TO  LIST  OUTPUT  FOR  SEGMENT (  1) 


in  1 

sTrainei 


Y 

N 

N 


DO  YOU  WANT  TO  LIST  OUTPUT  FOR  SEGMENT(  2) 

20  YOU  WANT  TO  LIST  PREBUCKLING  RESULTANTS  AND  RING  FORCES? 
DO  YOU  WANT  TO  LIST  FORCES  IN  THE  DISCRETE  RINGS,  IF  ANY? 
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2.  ACESNID  INPUT  DATA  FOR  ALL  CASES 

The  input  data  for  ACESNID  must  encompass  all  waves  N 
used  for  the  B0S0R4  model.  The  waves  are  specified  by 
NSTART-NFINISH. 


12  0  0  /NbRDER,NFENDS,NFCENT,NFCMPT,NOI 
S84E-5  5.  833E4  0.005  JRHOFL.  VSOUND.jhlR 
11110  111  /OUTPUT  FLAIjS 

1  if/ SETS  JSGBEG, JPTBEC, JSGEND, JPTEND  L 
1  2  27  /SEFS  JSGBE(i,JPTBE(i,  JSGEND  .JPTEND  CENTRA: 
115  /SOURCES  JSGB£G,JPTBEG, JSGEND. JPTEND  LEFT 
1  2  27  /SOURCES  JSGBE& , JPTBE6 , JSGEND , JPTEND  CYL 


CYLINDER 


3.  PIFLASH  INPUT  DATA  FOR  ALL  CASES 


BOSSY100-103  are  the  shell  mode  files  produced  from  the 
BOSOR4  runs.  Note  the  deleted  shell  modes  that  were  purely 
torsional  for  N  *  0. 


1  1  /NUMBER  NTORSN  NPTM  NSYM  NSYMP 

2  /(trWETSEG(K) , NUSESEG(K)  ,K=1,2 
/GRAVITY 
/SMF ,N*0 
7  SMF  11*1 
7 SMF ! Nr  * 


ivi 


?°8 


23  26  27  28  29  30 
19  20  21  22  23  24  23 

19  20  21  22  23  24  23 


.  -  7  ■  -  9  10  11  U  13  14  15  16  17  18 

VXA  YJA  7  8,9,10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25 
?6  27  28  29  30  /  (JUSE(J,K),  J*»JUSE(K)J 


4.  SAP  IV  INPUT  DATA  FOR  D  *  4.2  IN.  CASE 


The  following  data  file  generated  the  internal  structure 
model.  Note  that  only  half  the  doubly  overhung  beam  is 
modeled  as  indicated  by  the  beam  lengths  at  the  node,  i.e., 
node  21  ~  44.81250. 


■4.2  IN 
NRIGID 

.  0  / GEOMETRY  AND  DOF 


HI  756944 

U:°sim6 


9. QE*06 
9, 0E+06 

.4  0.0 


m 


8:8  ?  8:8 
0.  0  0.  0  0  0.  0 
o.  o  o.  o  i  o.  o 
0.  0  0.  0  0  0.  0 
MTENTS  TOTAL 


). 0/BEAM  MOD, POTS, DENS. 
5.0  /BRACKET  W 
35 /BEAM  AREA  AND  MOMENT 


1  8.4  0.00.0  1,02.8  12. 35/BEAM  AREA  AND  MOMENTS  OF  I 

2  6.3  O.O  O.O  l.  0  20.84  .  5 25 /BRACKET  m  * 

0.0  0.0  0.0  0.0 

o.  o  g.  o  g.  o  o.  g 

o. g  o.g,o. o  o. o  /  muult  c-x,y,z 
112  21  220  oooooo 
2  2  3  21  2  20000000 

l  l  ^2  21  2  g  0  0  0  0  000000  000001  0/SET  NODE  RELEASE  AT  12 
20  20  2i;,l  1  1  0  0  0  0  0  0  0 

5  g  g.  g  g.  o  o.  o  o.o  o.  o  o.  o  /conc  masses 
g. g  g. g  g. g  o.o  /element  load  multipliers 
9  9  9..2-.2_°:.9o  o.o  /ifpr.ifss,. ..  ,nfo, shift 


n  oAoVoV< 


NODE , f 
.(DUMlft 


TITLES 


5.  PICRUST  INPUT  DATA  FOR  D  3  4.2  IN.  CASE 


The  circumferential  location  of  the  substructure 
attachment  point  is  specified  by  ANDEG  which  is  the  angle 
ccw  from  the  shell  global  Z-axis.  The  meridional  location 
is  specified  by  the  node  21  at  frame  18.  The  angle  DEGROT 
specifies  the  orientation  of  the  substructure  local  z-axis 
to  the  shell  Z-axis.  This  is  indicated  by  in  Figure  5. 


1111011100 

0111110101  /OUTPUT  FLAGS  1-20 
500  /NJUSE,NHWSOB,NHWBAR 
1  2  21  180. 0/nIPSUB,LB0SEG,LB0SPT, ANDEG 
180.0  /  DEGROT 

12345/  (JUSE(J),J=1,NJUSE) 


6.  USLOB  INPUT  DATA 


In  Phase  III  of  ELSHOK  the  USLOB  code  performs  the 
integration  of  the  incident  pressures  using  a  Runge-Kutta 
scheme.  NCHARGE  specifies  which  type  of  UNDEX  is  to  be 
represented.  For  the  tapered  charge  nine  discrete  pressure¬ 
time  points  are  used  as  input.  The  numbers  following  the 
charge  weight  of  352  lb.  are  the  constants  used  in  equations 
(1)  and  (2)  that  empirically  describe  the  explosive  HBX-1 
(order-Kl,  K2,  K3,  K4) .  Note  the  0.0  for  surface  cutoff 
effects.  For  the  empty  shells  NSUBS  must  be  set  to  0.  The 
integration  is  performed  to  compute  velocities  at  all  the 
nodes  that  make  up  the  shell  and  substructure  but  only  a 
couple  are  requested  as  printed  output  in  NPTSH  and  NPTSUB 
as  indicators  that  the  program  ran  successfully.  Nodes  1/01 
and  2/21  are  checked  on  the  shell  and  node  1  and  21  are 
checked  on  the  substructure.  The  value  NQUAN  =  3  requests 
the  z  direction  velocities  as  output. 


a  / N SHAPE 

0.0  240. 0  .0002  320.0  .00053  350.0  .00243  305.0  .00294  160.0 
.00340  115.0  .0044  70.0  .00540  48.0  .00565  0. 0/DISCRETE  POINT  INPUTS 
0.0  /DECAY 

ioi8lOi010  /OUTPUT  FLAGS  1-20 


/NPTSHL 

1  2  21  /LBOSEG .LBOSPT 
/NPTSUB 

3  21  3  /NSTAT, NQUAN 


V..<u  «.  '■ 


•  -  "  .  •  .  - .  "j.  •  ,> 


n 


1 


.  t.-i  ijc  ».«•*,  f 


7.  POSLOB  INPUT  DATA 

This  data  file  generates  the  velocity  profiles  that  are 
plotted  on  Tektronix  devices.  The  time  step  DELT  and  the 
interval  of  integration  specified  by  NTIME  and  NSKIP  must  be 
the  same  as  for  USLOB.  This  input  requests  velocity 
profiles  be  expressed  in  terms  of  the  shell  global 
coordinate  system.  Shell  output  is  requested  by  segment 
node  and  angle.  Each  substructure  point  is  specified  by 
node,  number  of  degrees  of  freedom  at  the  node  for  which 
velocity  profiles  are  to  be  produced  and  which  degree  of 
freedom  it  is  produced  for.  The  tip  is  node  4,  one  degree 
of  freedom  in  the  local  z  direction. 

331  2  1  l  1 /NTIM , NSKIP .NSUBSjNTEK. HOARD 

VELOCITY^!)/ SEC* 

16/NPTSHL-W 

WivL 
n™iTio  uo.o 
I-f?lTo  1.0.0 
m  I'”.  I..,. 
mi'Sto  i.o.o 
rtmi.o  uo.o 
rtsvsio  uo.o 
rtrv{2a.o  i.«.o 

?-«V«.o  1.0.0 
vn  i'sL  uo.o 

V-ATH  2/01 

2  1  1  i.o.o 

Vfi  i  f  5.o 
rtiHiH15o.o 
rtivi  to.. 
raHi2iHo.o 

if 5U/LBoIeI , LBoIp?TnSPHC , NANO .ANGDEG 
16/NPTSHL-O 

m  i'8io 

rn  i'840  wo.. 

m  i'B’o  wo.. 


nr 


onnoa 


APPENDIX  B 


PROGRAMS  USED  TO  PREPARE  AMD  PLOT  VELOCITY  PROFILES 


PROGRAM  PLOTCONV 


THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  READ  THE  PUNCH- CARD 
FILE  FROM  PUSLOB  AND  CONVERT  IT  TO  A  FORMAT  WHICH  CAN 
BE  ACCEPTED  BY  EASYPLOT  OR  DISSPLA  ON  IBM  3033. 


DIMENSION  VELf 166.32) 
CHARACTER* 10  LaBeL 
character*4o  LTITLE 
COMMON^SOLVE/DTRECS 

DO  10~K=1.NCURV 
READ(  11,949}  LTITLE^ 


VELO( 166,2),  LABEL( 32) ,  TAG(32),  LTITLE( 32) 


>,R8fg  iDIBECS 


998 

997 


30 


wfej/jsi^REcs) 

F0RMATIA40 , I5.1PE11.  4) 

FORMAT'  ■  ’"’'-’fA  ' 

FORMAT 
DO  20 
T=0.  0 
LL=L*20 

DO  30  M=1,NRECS 
VELOfM.  lpT 
DTRECS 1=1000.  0*DTRECS 
T=T+DTRECS1 

My)=vEL(M’L)/u-° 


DO  40  KK=1,NRECS 

continue’”6* 


40  _ 

WRITEf  LL, .  . 

WRITEfLL  99 
20  CONTINUE 
996  FORMATf 2F15. 5) 

995  FORMATf A40 , 15 ' 1PE11.  4) 
994  FORMATf A10  Fl6. 5)  ' 


LTITLE 

LABEL( 


L*)!tAG(L) 


NRECS, DTRECS 


PROGRAM  FISH1  IS  THE  FORTRAN  CODE  REQUIRED  TO  PLOT 


VELOCITY  PROFILES  USING  DISSPLA 
REAL*4  VEL(I66),TIM(166),TIM2(I 
DO  10  1*1, 166 

DO  20  J=l, 166  „  ,  , 
READfal , *5TIM2(J)  ,VEL2(  J) 


166) ,VEL2( 166) 


20  CONTINUE 
CALL  COMPRS 


CALL  PAGE(8.  0,8.  0) „ 

CALL  AREA2D( 6.  , 5.  5) 

CALL  XNAMEf .TIME  CMSEC)$ '  1001 . 

CALL  YNAME{ ’VELOCITY  ( FT/SECT$ ! , 100 ) 

CALL  HEADIN( 'VELOCITY  PROFILE  $',100,2.  ,1) 


CALL  FRAME 
CALL  INTAXS 
CALL  YAXANGfO.) 
CALL  YTICKSt  5 ) 
CALL  XTICKS( 5) 


CALL  GRAF(0.  ,2.  0.10.  .-10.  ,10.  ,35.  )  „ , 

CALL  MESSAG('nODE  2/il  PORT  V-ATH  C$ ' ,23 , 1. ,5. 7 ) 
CALL  LINESP  (1. 75)  „ ,  v 

CALL  LINES  {'EMPTY  SHELLS' , IPAK.l) 

CALL  LINES  ('LOADED  SHELL$ ' , IPAK , 2 ) 

XW=XLEGND  (IPAK, 2) 


CALL  LINESP 
CALL  LINES  • 
CALL  LINES 
XW=XLEGND  ( 
YW=YLEGND  ( 
CALL  MYLGEN 
CALL  LEGLIN 
CALL  POLY3 
CALL  DASH 
CALL  CURVEf 
CALL  RESET t 
CALL  CURVE(' 
CALL  DOT  , 
CALL  GRID( 1 
CALL  RESETt 


$',100) 


CALL  RESElf ' DOT  )  „  ,  „ x 
CALL  LEGEND  (IPAK, 2, 4. 0,4. 5) 
CALL  ENDPLiO) 


CALL  DONEP 


AO-A167  893 


UNCLASSIFIED 


A  PARANETRIC  STUDV  OF  ELASTIC  RESPONSE  OF 
SU8HARINE-INSTALLED  EQUIPMENT  S.  .  <U>  NAVAL  POSTGRADUATE 
SCHOOL  HONTEREV  CA  G  F  DECONTO  HAR  86 

F/G  12/19. 1 


;r.  'r.dy.v.'.VCT \»  -.vr.i 


APPENDIX  C 

PROGRAM  USED  TO  ESTIMATE  AXIAL  FORCES  IN  THE 
SUPPORT  BRACKETS 


C 

c 

C 


PROGRAM  FORCES 
ESTIMATES  FORCES  DEVELOPED  IN  THE  SUPPORT  BRACKETS 


|AL*4  XL j, EMOD ! H, AREA ,VT0,VTllVT2f VBO , V& 1 , $BZ , DT | DB ) DEF 


REAL. 

REAL  ,7  AU  1  bi  JWV  I  U««*V 

rea^*4^xf6rc  ,  TtME ,  D 
EM^D=! 3E*08 

‘  fi 


DO  l8Ei=|,16S 

IKKiy  [ 3(5:81:8:1:1] 

10  CONTINUt 


DO  20  J=l, 165 ,2 
II*J 


JJ=II* 


VBC 

VB] 


0=VT 

rii,2) 

1=VT] 

JJ  2 

2=VTi 

KK  2 

0=VBl 

II  2 

l=VBl 

JJ  2 

2*VBi 

KK  2 

=  (  VTl 

3*4.  *V 

DB=(  VB0*4. 


*VT2)*H*4. 0 
__  w_.  .  „  __*VB2)*H*4.  0 

DEF=(DT-DB) 

?HR^iilS4?L)*EM0D*AREA 

$M5> 


20 


XDE^ 

88Nl5LL= 

WRITE56  * 
WRITE(_6  ,* 
25  CONTINUE  J.v 
WRITE( 83 . 


30 


fdJP’IH'1*2) 

nAAAUI  J  _  MSEC 

DO  38  E» 1.1*3. 2 
WRITE(83,*)  (jt(L,NN),NN 
CONTINUE 


STOP 

END 


TAPERED  CHARGE 


MSEC 

5. 9999999E-02 
0.  1800000 
0.  3000000 
0.  4200000 
0.  5400000 
0.  6600000 
0.  7800000 
0.  9000000 
‘ , 020000 
,  140000 
, 260000 
rJ  380000 
[.  500000 
,  620000 
740000 
360000 
980000 
100000 
220000 
340000 


2^0312 
-4871. 925 


♦  871.  ... 
-19920. 20 

-5310. 641 
5320. 271 
4192. 220 
7970. 726 

km 

-  p| 

24680. 4* 
8301.49 

m.\\i 


96 


LBF 


=  1,2) 


CONVENTIONAL  CHARGE 


MSEC 

5. 9999999E-02 
0. 1800000 
0.  3000000 


0. 4200000 
5400000 


0.  5< 

0. 6600000 
0.  7800000 
0. 9000000 

\m 

1.  380000 
‘  500000 
620000 
740000 
, 860000 
980000 
“3000 
3000 
3000 


62.  6$  _ 
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